Elementary excitations in dipolar spin-1 Bose— Einstein condensates 
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We have numerically solved the low-energy excitation spectra of ferromagnetic Bose-Einstein 
condensates subject to dipolar interparticle interactions. The system is assumed to be harmonically 
confined by purely optical means, thereby maintaining the spin degree of freedom of the condensate 
order parameter. Using a zero-temperature spin-1 model, we solve the Bogoliubov excitations for 
different spin textures, including a spin- vortex state in the absence of external magnetic fields and a 
rapidly rotating polarized spin texture in a finite homogeneous field. In particular, we consider the 
effect of dipolar interactions on excitations characteristic of ferromagnetic condensates. The energies 
of spin waves and magnetic quadrupole modes are found to increase rapidly with the dipolar coupling 
strength, whereas the energies of density oscillations change only slightly. 

PACS numbers: 03.75.Kk, 03.75.Mn, 67.85. Fg, 67.85.De 



I. INTRODUCTION 

Energetically low-lying collective excitations are cen- 
tral for understanding systems that exhibit superfluidity. 
For example, sound waves are the lowest-energy excita- 
tions in superfluid '^He, which explains why such a sys- 
tem may sustain frictionless flow below a critical veloc- 
ity, namely, the speed of sound. Similarly, the excitation 
spectrum of a superconductor has a gap in the neighbor- 
hood of the Fermi energy, allowing electric currents to 
flow without resistance. 

The study of collective modes in dilute atomic Bose- 
Einstein condensates (BECs) has had a key role in prob- 
ing the physics of these versatile quantum systems. Con- 
densates confined in strong magnetic traps are typi- 
cally described by scalar-valued order-parameter fields, 
and their low-energy excitations are density oscillations. 
The possibility of trapping BECs optically has liber- 
ated the atomic spin degree of freedom, providing an op- 
portunity for experimentalists to investigate condensates 
with spinor fields il!]. In the case of spin-1 alkali-metal 
BECs, such as ultracold vapors of ^^Rb in its ground- 
state hyperfine manifold, the order parameter is a three- 
component vector [2, 3] . The spin-1 gas of ^^Rb is magne- 
tized due to ferromagnetic interatomic spin coupling [j], 
and the related spin dynamics have been studied exper- 
imentally [a, Q . The lowest-lying collective excitations 
constitute a diverse set of density oscillations, spin waves, 
and quadrupolar spin fluctuations. 

Recently, long-range anisotropic dipolar forces between 
the condensed bosons have attracted considerable inter- 
est. It has been theoretically predicted that dipolar in- 
teractions can give rise to interesting magnetic structures 
such as spin vortices [^-Q and helices Ifljllll- Also, dipo- 
lar forces have been suggested to be involved in the dis- 
integration of a helical spin texture into small magnetic 
domains in a condensate of spin-1 *^Rb , 12;] . In a similar 
system, dipolar interactions have been shown to induce 
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magnetoroton softening in the excitation spectrum, link- 
ing the study to the physics of ^He [ij] . In this article, we 
investigate how dipolar interactions affect the elementary 
excitation spectrum of a harmonically confined pancake- 
shaped spin-1 BEC. In particular, we solve the Bogoli- 
ubov spectra of a spin- vortex and a nearly spin-polarized 
(flare) state in the absence of polarizing external fields. 
In addition, we calculate the excitations of a dipolar BEC 
with the atomic spins polarized perpendicular to an ex- 
ternal field. In a scalar dipolar system, collective exci- 
tations have been studied both theoretically fl4l-[l6| and 
experimentally [17!|. 

The stationary states studied here using a spin-1 model 
have analogous counterparts in a semiclassical framework 
that describes condensates of particles with large dipole 
moments d, [lal- Thus, the presented results may shed 
light on the behavior of strongly dipolar systems such as 
the spin-3 condensate of ^^Cr ^9i, t20,] which has atomic 
magnetic moments of 6 /xb and which has been produced 
by purely optical means |21| . Other strongly magnetic 
atomic gases that have been successfully cooled to ul- 
tralow temperatures include vapors of Tm [22[, Er j23| . 
and Dy |24| . which have atomic magnetic moments of 
4/iB, 7/xb, and 10 /iB, respectively. 



II. THEORY 

In this work, we study spin-1 BECs using a zero- 
temperature mean-field treatment, thereby neglecting 
possible effects due to the presence of a thermal com- 
ponent. The state of the condensate is described by 
a spinor wave function ^(r,t) = (?/'i,?/'0: "0-1) in the 
Zeeman basis. Its time dependence is governed by the 
Gross-Pitaevskii equation (CPE) 

a 



::h^- V 



Here, h^^ — — fi,^V^/2m + T4rap(r) — pSz + qS^ is the Bogoliubov equations of the general form 
single-particle Hamiltonian with the optical trapping po- 
tential Vtrap = m [ujf.{x^ + y2) + u}^z^] /2. The last two f A B 
terms of h^^ account for the linear and quadratic Zee- \—B* —A" 
man terms with the external homogeneous magnetic field 
pointing along the z axis. The constant p depends lin- 
early and q quadratically on the strength of the applied 
field. In Eq. ^, n(r) = *t^ ^ EL-i V'fcV'fe is the parti- 
cle density and Ma(r) = 'ii'^Sa'^, a £ {x, y, z}, stands for 
the ath component of the local magnetization with {Sa} 
denoting the standard spin-1 operators. Furthermore, 
la (r) is a spin-dependent effective potential arising from 
nonlocal dipolar interactions (see Appendix). The cou- 
pling constants determining the strengths of the density- 
density, spin-spin, and dipole-dipole interactions are 
given by g^^ = Anfi^iaQ + 2a2)/3m, gP = 'iTTh^{a2 - 
ao)/3m, and g^ = no^^gp/An, respectively. Here, /ig 
denotes the permittivity of the vacuum, fiB is the Bohr 
magneton, and gp is the Lande factor (gp = ^1/2 for 
s'^Rb with F = 1). 

From now on, we will assume that the components of 
the order parameter have the same Gaussian profile along 
the z axis, *(r) = 'i'{x,y)A^exp{-z'^/2<7^) with Aa- = 



(4) 



CT^^'^TT^^'^, implying that the direction of magnetization 
does not depend on z. This approximation should be 
accurate in the limit of tight axial trapping, lUz ^ uir- In 
the noninteracting limit, the ansatz produces the ground 
state if cr = Oz = y^h/mtUz- Thus, in the case of a 
pancake-shaped BEG, we may reduce the GPE to the 
effectively two-dimensional form 



ihdt'^{x, y,t) ^ \h + gnU + ffs ^ M^Sa 

a. 

+ gd^/:^a}*(x,j/,i), (2) 



where the effective potentials are now evaluated at z = 0, 
h = -h^ {dl + dl) /2m -I- mujl (x2 -|- y^) /2 _ pSz + qSl 

and gi — gf^/\/2na, i € {n, s, d}. The dipolar integrals 
{/q (x, y)} include the width a of the Gaussian profile as 
a parameter. The precise forms of these functions are 
derived in the Appendix. 

Stationary states are obtained with the ansatz 
^{x, y, t) = ^{x, y) exp {—i^t/h), where /x is the chemical 
potential determined by normalizing the wave function to 
the total particle number N ^ J dr'^^^. In the presence 
of strong axial confinement, excitations of the conden- 
sate are frozen in the axial direction. Small-amplitude 
oscillations about a stationary state '^{x, y) are found by 
assuming the wave function has the form 

*(a;, y, t) = e"*^*/'^ h{x, y) + u(x, y)e-"^' + w*{x, y)e"^'* 

(3) ■ 
where the vectors u and v are called quasiparticle am- 
plitudes. Substituting this trial into Eq. ^ leads to the 



Here, hw denotes the excitation energy of the mode 
{u,v}, and A and B are linear operators whose detailed 
forms are given in the Appendix. 

An external magnetic field B exerts a torque M x B 
on a magnetic moment M, which causes the moment to 
precess about the direction of the magnetic field with 
the Larmor frequency wl. Thus, in the case of a ho- 
mogeneous axial field BzSz, the magnetization in a sta- 
tionary state is polarized along the z axis. Yet, at 
least in the absence of dipolar interactions, we can allow 
for more general magnetization textures by considering 
states that are stationary in a spin reference frame that 
rotates about the magnetic field. Such states are of the 
form "^p = exp(ipSzt/h)'i!o, where p/h is the precession 
frequency about Sz and ^q is a stationary solution of 
Eq. ^ with p = 0. This simple picture holds when the 
dipolar interactions are absent, but when gj ^ 0, the 
dipolar potential is time dependent in the rotating frame 
and the analysis becomes more involved. However, for 
systems considered in this work, the Larmor frequency 
in a typical magnetic field of strength 100 mG is in the 
MHz regime, which is much faster than other relevant 
dynamics of the condensate, such as spin and density os- 
cillations whose characteristic frequencies are in the kHz 
regime. Therefore, the system can be viewed as being 
subject to an effective dipolar potential obtained by time- 
averaging the instantaneous dipole-dipole potential over 
one Larmor cycle [lOj . The time-averaged potentials are 
given in the Appendix. 

In the numerics, we discretize the effectively two- 
dimensional equations on a uniform grid. The Bogoli- 
ubov matrix in Eq. ^ is typically full when g^ > due to 
nonlocality of the dipolar potentials, and explicit diago- 
nalization of the large eigensystem becomes cumbersome. 
Instead, we use the implicitly restarted Arnoldi iteration 
implemented in the ARPAGK library for solving only the 
relevant part of the spectrum, namely, the energetically 
low- lying eigenmodes. The method relies on repetitive 
operation by the system matrix without the need to store 
it in memory. The two-dimensional Fourier transforms 
used in evaluating the dipolar potentials are implemented 
efficiently with a fast Fourier transform (FFT). 

In the following, we choose a^ = y^h/mujr and htOr for 
the units of length and energy, respectively. We set the 
dimensionless coupling constant to g'^ = g^^ N / HuJraf. = 
500, which corresponds roughly to 10"* ^^Rb atoms in 
a trap with uj^j^-n — 100 Hz. Moreover, we choose 
a = a/2/5, yielding ~g^ = ~g^^ I\r2^a k. 700. The conden- 
sate is assumed to be ferromagnetic with g^jg-a = —0.01. 
Qualitatively, the results we present are not very sensi- 
tive to the values of g^, gs, or a. Instead of fixing the 
dipolar coupling strength to some predetermined value 
(for ®^Rb, gd/gn = 4.2 x 10~^), we consider a range of 




FIG. 1. (Color online) Schematics of the three spin textures 
considered in this work: (a) spin vortex, (b) flare, and (c) 
spin-polarized state. The spin vortex and flare are studied in 
the absence of magnetic fields, whereas the spin in the polar- 
ized state is precessing rapidly with the Larmor frequency wl 
about a homogeneous field perpendicular to the plane of the 
figure. 




FIG. 2. (Color online) Total energies of the spin-polarized 
state in an external axial magnetic field (solid line) and the 
flare state (dashed line) as a function of the dipolar coupling 
strength g^. The energies are given with respect to the total 
energy of the spin-vortex state shown in the inset and by 
the dotted line in the main flgure. The Mermin-Ho vortex 
(dash-dotted line) is an excited state for all parameter values 
considered. 



interaction strengths within the entire stability window 
< gd/gn < 0.24 [i,[23. Although alkali- metal conden- 
sates are typically subject to weak dipolar forces, it might 
be feasible to tune the ratio gd/gn experimentally by uti- 
lizing optical Feshbach resonances [2^ [231 . For the exter- 
nal magnetic field, we choose B = B^Sz with Bz = 100 
mG. In such a field, the Larmor frequency of a ^'^Rb atom 
is ojl = gpfJ-sBz/fi ~ 27r X 70 kHz, and the quadratic Zee- 
man shift has a magnitude of q/h sa 27r x 0.70 Hz |28j . 



III. RESULTS 

The homogeneous spin-polarized spin-1 system has a 
rich variety of excitations 0, Q: The collective den- 
sity oscillations obey the typical Bogoliubov dispersion 
relation e(fc) = y/e^ik) + 2gnneo{k), where eo(fc) = 



h k /2m. The long-wavelength density modes have a lin- 
ear phonon-like spectrum and the short-wavelength exci- 
tations approach the quadratic free-particle case. The 
spin waves have a free-particle-like spectrum e{k) = 
£o(fc). The magnetic quadrupole oscillations are pe- 
riodic modulations in the quadrupole tensor Qa/B = 
2Sap/3 — (SqSg + SpSa)/{2n) that describe local spin 
fluctuations [29|. They have a quadratic dispersion with 
an energy gap, eik) — SQik) + 2\gs\n. Here, we consider 
these excitations in an inhomogeneous system in the pres- 
ence of dipolar interparticle forces. 

We have solved the low-lying excitations in a ferro- 
magnetic dipolar spin-1 BEG for the three types of sta- 
tionary states illustrated in Fig. [1] First, we consider 
two stationary states in the absence of an external po- 
larizing field;namely, the spin vortex [Fig. [TJa)]and the 
flare [Fig. [Hb)], studied previously in Refs. 041 [ll]. 
We choose the spin-vortex order-parameter field to be 
of the form ipkix, y) = tpk{r) ex.p{—ik(p) . In such a state, 
the transverse magnetization rotates through an angle 
of 27r about the unmagnetized vortex core. Second, we 
consider the ground state with vanishing axial magne- 
tization in a homogeneous axial magnetic field. Due to 
the rapid Larmor precession of the transverse magnetiza- 
tion, the condensed particles feel a time-averaged dipolar 
potential, which results in a transversally spin-polarized 
texture [Fig. [T{c)]. The total energies of these station- 
ary states as a function of the dipolar coupling strength 
5d are depicted in Fig. [51 The inset shows the energy 
of the spin-vortex state, which is the ground state for 
gd/gn > 0.024. For comparison, we also show the energy 
of the Mermin-Ho vortex that carries both mass and spin 
currents. 

The spin vortex and the spin-polarized state are cylin- 
drically symmetric, and therefore their excitations can be 
written in the form 

u,(r,v,)=e<''-^^-)'^u,(r), (5) 

v,(r,^)=e<^+^^^)^v„(r). (6) 

We have chosen the spin winding s = 1 for the spin vor- 
tex and s = for the spin-polarized state. The quan- 
tum number k € Z determines the phase winding of the 
excitation. We are interested in low-energy excitations 
which typically have small values of n, such as breath- 
ing {k = 0), dipole {k = ±1), and quadrupole {k — ±2) 
modes. For fixed k, there are several low-energy exci- 
tations depending on the nature of the mode (density 
wave, spin wave, or magnetic quadrupole wave) and the 
number of radial nodes in the components of u^ (r) and 

Figure [3| shows the dependence of the excitation ener- 
gies on the dipolar coupling strength gj- In the limit of 
vanishing ^d, the low- lying parts of the spectra host myr- 
iads of excitations, of the order of a hundred in the range 
£ G [0, 3 huir] for each state. However, in the limit of 
strong dipolar interactions, only a few excitations remain 
within this energy window. For clarity, we have presented 
the energies only for the most interesting modes. 
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FIG. 3. (Color online) Excitation energies of selected Bogoliubov modes as a function of the dipolar coupling strength gd- The 
graphs refer to the (a) spin vortex, (b) flare, and (c) spin-polarized state illustrated in Fig. [l] In the limit gd/gn — >■ 0, the 
three highest-energy modes in each figure are the density breathing (k = 0), quadrupole (k = ±2), and dipole (k = ±1) modes, 
in order of descending energy. The three lower-lying modes are the lowest magnetic quadrupole oscillation (k = 0) and spin 
waves (k = ±2 and k = ±1). The energies of the magnetic quadrupole mode and the spin waves increase rapidly with dipolar 
coupling strength, whereas the energies of the density oscillations are less affected. The density and spin oscillations become 
mixed due to the dipolar forces. In the spin-vortex case, the density and spin oscillations of similar symmetry exhibit avoided 
crossings. 



In Fig. ^a,) , the three highest-lying modes in the limit 
gd/gn -^ are the density breathing, quadrupole, and 
dipole modes with energies (2, \/2, l) hur and quantum 
numbers k = 0,±2,±1, respectively. The excitation en- 
ergy of the dipole mode is determined solely by the ex- 
ternal trapping potential, and hence it is independent of 
internal interactions, such as dipolar forces. In general, 
excitations that induce strong density fluctuations are 
relatively robust against changes in the dipolar coupling 
strength. The three remaining modes shown, in order 
of decreasing energy, are the lowest magnetic quadrupole 
oscillation (k = 0), the quadrupole spin wave (k = i2), 
and the core-localized mode (k = ±1). The energies of 
these excitations increase rapidly with increasing gd/gn- 
The core-localized excitation corresponds to an imagi- 
nary eigenvalue for gd/gn ^ 0.01, which signifies the dy- 
namical instability of the spin vortex. 

The spin-vortex state ipk{x,y) = ■0/c('') exp(— zfc(/j) is 
invariant under the discrete symmetry ipk -^ ''P-k- ^^~ 
citations without this invariance are doubly degenerate. 
In Fig. inta), these include all but the density breath- 
ing mode and the lowest magnetic quadrupole oscillation. 
For the Mermin~Ho vortex 4'k{x, y) = ipkif) exp[— z(A: + 
l)ip], this symmetry is broken, and hence the spectrum is 
mostly nondegenerate. The two core-localized modes of 
the Mermin-Ho vortex have negative excitation energies 
for weak dipolar interactions, gd/gn ^ 0.1, which implies 
that the Mermin-Ho vortex is locally energetically un- 
stable. 

The flare state is only found within the interval < 
gd/gn ^ 0.07, beyond which we obtain a state with the 
same symmetry but hosting two spin vortices. Ener- 
gies of several excitations of the flare state are shown in 
Fig.|3Ib). The three highest- lying excitations in the limit 
gd/gn — > are the density breathing, quadrupole, and 



dipole modes. The two modes below them are the lowest 
magnetic quadrupole oscillation and the dipole spin wave. 
The dipolar interaction breaks the rotational symmetry 
of the flare state, which lifts the degeneracy of all but 
the density dipole mode. Previously, the frequency split- 
ting of the counter-rotating density quadrupole modes 
due to the presence of a mass vortex in a scalar sys- 
tem was utilized to measure the angular momentum of a 
BEC [30, |3l|. In the absence of external magnetic fields, 
the lifted degeneracy in Fig. ^ih) could provide a spec- 
troscopic means of measuring the strength of the dipolar 
forces. For gd/gn ^ 0.05, the lower branch of the dipole 
spin wave turns imaginary. 

In Fig. ISjc), the excitation energies for the transver- 
sally spin-polarized state in a homogeneous axial field 
are shown analogously to the previous cases. The three 
lowest-lying modes in the limit gd/gn ~> are the lowest 
magnetic quadrupole oscillation (k = 0) and two lowest 
spin waves {k = ±2 and k, = ±1). 

The avoided crossings of the excitation energies occur 
for density and spin waves of the same quantum number 
k; for example, such a crossing is visible in Fig. El^a) 
for the quadrupolar density and spin oscillations with 
K — ±2 at gd/gn ~ 0.06. A similar phenomenon has 
been previously studied in anisotropically trapped scalar 
BECs [22, ISJl • In our case, the perturbing Hamiltonian 
is the dipolar interaction term. The interaction couples 
the neighboring quantum numbers k in the spin-vortex 
case but does not have an azimuthal dependence in the 
rapidly precessing spin-polarized state. Therefore, the 
avoided crossings appear for the spin vortex in Fig.[31[a), 
as well as for the Mermin-Ho vortex, but are absent in 
the spin-polarized case. 

In the absence of dipolar interactions, it is usually 
straightforward to distinguish the density waves from 




FIG. 4. (Color online) Condensate density modulation An 
caused by a given excitation of the spin vortex compared to 
that of a simple translation Ano. In the limit gd/gn — >■ 0, the 
curves show the modulations due to the density breathing 
(solid), density quadrupole (dashed), density dipole (dotted), 
and the quadrupole spin (dash-dotted) modes. The roles of 
the quadrupole density and spin waves are interchanged in 
the neighborhood of the avoided crossing shown in Fig. [3{a). 



spin waves in the excitation spectra. The long-range 
anisotropic dipolar forces naturally couple spin and den- 
sity oscillations. In order to quantify the change in the 
density due to a given excitation, we define the time- 
averaged integrated density modulation jsj 



IV. DISCUSSION 

In summary, we have numerically studied the collec- 
tive excitations of a harmonically confined dipolar spin-1 
condensate. The Bogoliubov quasiparticle spectra were 
solved for three types of spin configurations often encoun- 
tered in ferromagnetic systems: the spin vortex, the fiare, 
and the spin-polarized texture precessing rapidly due to a 
homogeneous external magnetic field. Particular empha- 
sis was placed on investigating the effect of dipole-dipole 
interactions on the quasiparticle spectra: the excitation 
energies of spin waves and magnetic quadrupole oscil- 
lations were found to increase rapidly as a function of 
the dipolar coupling strength, whereas the energies of 
excitations related to collective density oscillations were 
essentially unaffected. 

As an extension of the present work, one could solve 
the excitations of dipolar pancake-shaped condensates 
in a wider energy range in order to study the effect of 
trapping on the roton-maxon spectrum. Previously, the 
emergence of the roton minimum has been investigated 
in the quasi-two-dimensional limit both in scalar [35] and 
spin-1 systems [13]. Also, taking the dipolar forces into 
account in studying the dynamical instability of the spin 
helix of Ref. [3a] mig ht be relevant due to their possible 
stabilizing effect |11| . 

It would be interesting to investigate the collective ex- 
citations at nonzero temperatures. However, based on 
earlier studies of scalar condensates [33, |3^ , we expect 
that finite-temperature effects do not alter the qualita- 
tive features of the excitation spectra at temperatures at 
which experiments with dilute BECs are typically carried 
out. 



An= /'dr|(*'^1''-*t^)| = -e /"dr|*^u + *'^v|. (7) 
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Appendix: Effective dipolar potentials 

Here, we briefly present the steps required for deriv- 
ing the effectively two-dimensional GPE, Eq. ([2]), and 
calculate the dipolar integrals {I^{x,y)} therein. Re- 
peated indices are to be summed over throughout the 
Appendix. All but the last term in Eq. ([2]) are obtained 
in a straightforward manner by substituting the ansatz 
*(r) = ^{x,y)G{z) with G{z) = A„exp{-z'^/2a^) and 
A„ = cr"i/27r-i/4 jj^^Q gq (p^ multiplying by G{z), and 
integrating over z. The dipolar integrals in Eq. ([1]) read 

/„(r) = f dv'D^p{v - Y')Mfi{v'), (A.l) 



with Dai3{r) = (Sapr"^ — Svarp) jr' . After the aforemen- 
tioned substitution procedure, the dipolar term becomes 
9Ala{x,y)Sa^{x,y), where 



il^^.v) = ^ fdr'D^pir-r')Mp{x\y')Giz') 

(A.2) 
The right-hand side is easily simplified with the convolu- 
tion theorem. The required Fourier transforms are given 
by Hi 



47r 



Arr 



By denoting M^ — J^jM^j}, we obtain 



^"^^>=if„Mk)(A.3) 



(A.4) 



i:i^,y) 



I-kgT 



\i^aSi\Kxii^y,i^z) 



xM^(fc,,fc,)e-('==^)'/^} 



(A.5) 



2 = 



Because the expression is evaluated at z = 0, the inverse 
Fourier transform over fc^ simply reduces to an integral. 
By evaluating the integrals and denoting 



KL. = ^/(fc.a) = K' 



KL 



^y jL2 •' ^ ' ' y^ 

An kl 



r|/(fc.a), 



Kz = Y -/(fcrCr), 

-^a^ = otherwise, 



(A.6) 

(A.7) 
(A.8) 

(A.9) 
(A.IO) 



where /(O = (27r)^/^^exp (fV2) erfc (^/^/2) and fc^ = 



1,2 



fc^ , we find 



/^(x,2/) =^-i{A';.^(fc„fc,)M^(fc,,fc,)}. (A.ll) 

Hence, within the assumption of the Gaussian axial pro- 
file, updating the dipolar potentials for a given state 
^(x, y) requires the evaluation of three two-dimensional 
Fourier and inverse Fourier transforms. By defining the 
linear operator L^pi-} = J-^JAT^^J"!-}}, Eq. ((XTT|) 

can be written simply as /^ = Cap{Mp}. 

In the presence of a homogeneous external magnetic 
field BzGzi the magnetization of the BEG precesses about 
the z axis with the Larmor frequency. The time-averaged 
effective dipolar potential is obtained from the instanta- 
neous form, Eq. (|A.11|) . by setting K',^y = and replacing 
K'^, and K'yy with {K'^^ + K'yy) /2 =_-KJ2. Thus, in 
the special case of a condensate polarized in the x direc- 
tion, the time-averaged effective dipolar potential takes 
the simple form Vd — —gdCzz{Mx}Sx/2. 

With the notation introduced above, the two- 
dimensional time-dependent GPE can be restated as 



+ gda0{'^^S0-i'}So)-i'. (A.12) 



Direct substitution of the trial ^(r,i) — 
e"'^*/'* [*(r) -Hue"*"* -|-v*e''"**] reveals that, in 
the absence of dipolar interactions, the linear operators 
A and B appearing in Eq. (j4]) can be written as 



yVn 



h + .g„ (*t^ + **^) 
+.9s (*^5„*^„ + ^„**t^„) , (A.13) 

(A.14) 



B„, = 5„**^-fg,^„**^^, 



T Q* 



When dipolar interactions are present, the operators A 
and B are augmented with the nonlocal terms 

id = .gd (^Cap{-^^Sp^}S^ + £^p{^^Sp ■ }5„4()V15) 
Bd = 5d>C„^{*^5^ • }5o*. (A. 16) 

Therefore, in the general case, the operators appearing 
in Eq. (g]) are given hy A = Ans + Ad and B = B^s + Bd- 
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